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Abstract 

We investigate the long-time limit of quantum localization of the kicked Rydberg atom. The 
kicked Rydberg atom is shown to possess in addition to the quantum localization time tl a second 
cross-over time t^) where quantum dynamics diverges from classical dynamics towards increased 
instability. The quantum localization is shown to vanish as either the strength of the kicks at 
fixed principal quantum number or the quantum number at fixed kick strength increases. The 
survival probability as a function of frequency in the transient localization regime tl < t < to 
characterized by highly irregular, fractal-like fluctuations. 

PACS numbers: 32.80.Rm, 05.45.Mt, 72.15.Rn, 05.45.Df 
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I. INTRODUCTION 



The quantum mechanics of classically chaotic few-degrees of freedom systems has become 
an intensively studied in the field of "quantum chaos" One key feature is quantum 

localization, i.e. the localization of the quantum wavefunction while the corresponding 
classical distribution shows diffusion . In periodically driven systems, this effect has 

primarily been studied in the kicked rotor, e.g. Q, S, 0, El, and in the Rydberg atom in a 



1 tne KicKed rotor, e.g. 



sinusoidal electric field, e.g. j8|, |9|, [lOl, [111, [l^, ll3[. For both systems, quantum localization 
is closely related to Anderson localization in transport in disordered systems One 
signal of this analogy are strong fluctuations in the quantum system for different observables 
when varying some external parameter. A third system ex per imentally and theoretically 
studied is the periodically kicked Rydberg atom (e.g. ^HHH 3, IS, 20, 21, 3, 3)) 
i.e. a hydrogen like atom prepared in a high-lying state and subjected to a sequence of ultra- 
short impulses. In a recent publication, we showed the existence of quantum localization 
in the positively kicked Rydberg atom i-e. the hydrogenic system with the initial 
state prepared in a highly elongated quasi-one dimensional state localized on one side of the 
nucleus and the periodic impulsive momentum transfer Ap > pushing the electron away 
from the nucleus. The classical phase space of this system is globally chaotic with all tori 
destroyed for arbitrarily small Ap, and the classical survival probability decays algebraically, 
Psur oc t''^, with a ~ 1.5 (see Fig. [T]and [3|)- contrast, by following the time-dependent 
system up to a few thousands of kicks we could show both the quantum suppression of 
classical ionization and the "freezing out" of the wavefunction, the hallmarks of quantum 
localization. Several issues remained open, however. They include 1) The origin of a slow, 
yet noticeable, decay of the localized states, 2) the dynamical role of high harmonics (up to 
infinity) present in the system, and 3) the properties of strong fluctuations present in the 
localization regime. 

In the following paper we address these issues. We identify two characteristic time-scales 
("break times") in this system. In addition to the localization time r/, where quantum and 
classical dynamics begin to differ from each other due to universal destructive interferences, 
there is a second break time td where localization is broken. Beyond td, a second cross-over 
occurs where the classical dynamics becomes more stable. This second cross-over in the 
kicked atom is related to the presence of non-classical photoionization. 
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In the next section, Sect.|TTl we describe the method used in our studies and in Sect. lIIII we 
show that quantum locahzation in the kicked Rydberg atom is transient. Strong fluctuations 
in the survival probability as a function of frequency for fixed times are studied in Sect. IIVI 
and in the last section (Sect. a summary is given. 



II. METHOD 



The Hamiltonian of the one- dimensional kicked Rydberg atom is (in atomic units) 



K 



H{t)=Ho-qApJ2^it + T/2-kT) , (1) 



k=l 



where Hq = ^ — Ms the hydrogen Hamiltonian and q the position of the electron. Ap and 
T are the strength and the period of the train of kicks, respectively. We will use the number 
K of kicks and time t = K xT interchangeably. The restriction to a ID model in the present 
context is necessary since the study of the long-time limit in 3D is currently computationally 
not feasible. Our previous studies for up to 10^ kicks have shown that the ID model can 
reproduce essential features of the 3D problem. For a more detailed discussion of the relation 



24| and references therein. In the 



between the ID model and the real 3D dynamics, see 
case Ap > 0, i.e. the kicks directed away from the nucleus, the classical phase space is 
void of stable islands and the effect of quantum localization has been shown to set in within 



a few hundreds of kicks 



24| . In the opposite case Ap < stable islands persist allowing 



for survival of Rydberg states both classically and quantum mechanically, referred to as 
stabilization j22| . 

The unidirectional kicks build up an average field F^'^ = —Ap/T. Hence, the time- 
periodic Hamiltonian (P) can be decomposed into the time-independent Stark Hamiltonian 



i^stark = ^0 + qF^-" (2) 
plus an infinite series of harmonics of equal strength. 



H{t) = i/stark + 2F'^^g COS (^27rm/T (^t - 



(3) 



In this paper, we will highlight the influence of the higher harmonics m > 1, their presence 
distinguishing our system from the Rydberg atom driven by a microwave field. For the 



"positively kicked" Rydberg atom with Ap > (i.e. < 0) the Stark Hamiltonian 

possesses a potential barrier with maximum at E^'^"^^^ = — a/2 | | and q^''^^^^'^'^ = a/ 1 / 1 F^^ \ 
resulting in a finite number of quasi-bound states and a continuum. To keep the Stark 
Hamiltonian invariant, i.e. the average field fixed, we vary u and Ap such as to keeping the 
average field fixed. 

To calculate the long-time evolution of the quantum system we represent the period-one 
time-evolution operator 

f/(T) = exp(-iiJo T/2) exp{iApq) exp{-iHo T/2) (4) 

in a basis | n) defined by Hq \ n) = En \ n) by means of the pseudo-spectral method 
Dirichlet boundary conditions are applied at g = and q = gmax- Solving the eigenvalue 
equation U{T) \ (pj) = exp{—iTSj) \ (pj) yields the time-dependent wavefunction in terms of 
Floquet states | as 

I i^{KT)) = U{T)^ I V^(O)) = expi-iKTS^) \ 0,) (5) 

j 

with dj = ((t>i I ijj{0))- A masking function in q is introduced to avoid spurious reflections at 
?max (see j2J]). For the low frequencies used in this paper, we apply the masking operator 
three times per period. The convergence of the wavefunction obtained is tested both by 
varying g^ax and by comparison with direct solutions of the time-dependent Schrodinger 



25| 



equation (i. e. without facilitating the Floquet states (f)j) as described in 

For later reference we introduce scaled units, denoted by the subscript 0. They leave the 
classical dynamics invariant and are defined hj Eq = E x nf, To = T/(27Tn^), uq = I/Tq, 
and Fq^ = F^^ x nj where rti is the principal quantum number (action) of the initial state 
[lol | . Note that the transition energy due to to absorption of one photon E^ = nf x hw = 
nf X 27r/T = vt^jui^ is not scaling invariant. 



III. TRANSIENT QUANTUM LOCALIZATION 

A. Survival probability and effective quantum number 

One measure to study quantum localization is the survival probability, defined as 
Psur(-^) = (V'(-^r)l-PboundlV^(-^r)) whcrc Pbound IS the projcctiou operator onto bound hy- 
drogenic states. The underlying picture is that transport along the energy axis towards 
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the continuum is considered to be the equivalent of conductance in disordered systems j3| ■ 
Correspondingly, suppression of energy absorption from the pulse sequence and, thus, sup- 
pression of ionization is identified as localization in a purely chaotic system or stabilization 
when the classical phase space is mixed. 
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FIG. 1: Classical P^^^ (dashed line) and quantum (solid lines) survival probability for the 
positively kicked Rydberg atom as a function of number K of kicks for |i^o^| = 0.005 (Apo ~ 0.021). 
The quantum data are shown for 11 frequencies z^o uniformly distributed between 1.45 and 1.45008, 
Hi = 50. (a): Long-time behavior, the dashed-dotted line indicates a fit to a power-law with 
exponent a = —1.5, Eq. and (b): short-time behavior. 
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FIG. 2: Mean principal quantum number (n), Eq. for the same parameters as Fig. Q {n)ci' 
dashed line; {n)qj^: solid lines, (a): Long-time behavior (of log(n)), the dashed-dotted line in- 
dicates a fit to a power-law, and (b): short-time behavior of (n). The arrow in (b) indicates an 
approximate localization time defined in Eq. ()16|). 

The classical P^^j.{K) is obtained by the classical trajectory Monte Carlo method (CTMC) 
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with a micro-canonical ensemble to represent the initial state |27l . |28| . The classical phase 
space (for Ap > 0) is fully chaotic without any tori left intact. We show an example for the 
survival probability in Fig. ^ The classical P^^^{K) for times larger than Kq kicks decays 
algebraically jlS^ . 

P:1{K) = (K/Ko)-- (6) 

with a ~ 1.5. Both the time Kq and the factor a are only weakly dependent on the 
parameters of the field (z/q and Apo)- By contrast, the quantum survival probability P^{K) 
displays a very different and much more complex behavior. One intriguing feature is an 
extreme sensitivity to the driving frequency uq, especially in the long-time limit (compare 
also to the results in e. g. 0,0, [ll|). The resulting fluctuations in the kicked atom will be 
studied below in more detail (see Sect. II V|). 

For characterization of quantum (de)localization it is useful to introduce another ob- 
servable that describes the bound portion of the wave packet. Inspired by the relation 
n = l/\/—2E for bound states, we introduce a mean quantum number as 

{n{K)) = (^(irT)lPbound ^4= Pbound|^(i^T))/P,,,(ir) . (7) 

In this equation, the expectation value is calculated for the portion of the wavefunction 
residing in the bound subspace. Equation (jZj) characterizes the specific position of that part 
of the wavefunction that remains bound, in particular that remains localized or stabilized 
near the initial state (or initial torus). Indeed, (n) displays a characteristically different 
behavior classically {n)ci and quantum mechanically (n)qm, both at short times (Fig. ^p) 
and long times (Fig. |2K). The sharp drop of (n)qm from 50 to 1 at ~ 10^ will be analyzed 
in Sect. IIII (^1 The corresponding non-normalized distributions P{n,t) whose mean values 
are given by (n) are depicted in Fig. IHl The quantum P^™(n, t) remains well localized close 
to the initial state while the time evolution of P^^{n,t) features a rapid increase of width in 
the energy distribution and a shift towards increasingly lower energies (i.e. smaller n). 

Note, that classically n does not possess a lower bound n = 1, as the quantum distribution 
does. The fact that the energy distribution can grow to negative values without bound is 
at the root cause for the algebraic decay. The probability for crossing the E = Q line (i.e. 
the n = oo line) in the next time step and thus for reducing P,^j. decreases as the mean 
distance to the = line as measured by {n{K))c\ increases with t. We now relate the 
power-law decay Eq. ^ to the motion of P'^^{n,t) away from the threshold. The classical 
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FIG. 3: (Color online) Time-dependent spectral distribution P{n,t), both classically (a) and 
quantum mechanically (b) for the same parameters as Fig. [21 (data shown for uq = 1.45). 

survival probability for large times is well fitted by Eq. ^ with Kq ^ 1000. The scaled 
average quantum-number is also well described by a power-law, 



(no(K) )ci = KK/Ko) 



-0.5 



(8) 



where the power-law decay is assumed to start at the same time as that of the survival 
probability and 6 ~ 1.8. The differential form of Eq. (jH)), 



dK 



j^^ surV-"- 1 -I 



(9) 



can now be expressed by means of Eq. (jH)) as follows, 



2 DCl 
sur 



(10) 



dt Kq 

with d = a/b"^. For some cases with stronger average fields (Fq^^ = 0.02 and 0.05) shown 
in Figs, m and El (where Kq ^ 60 and 9, respectively) we find that the power-law for (no)ci, 
Eq. (|H1), holds with approximately the same constant b. Hence, Eq. fjlOj) is still applicable 
and the constant d has essentially the same value. 
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FIG. 4: (n) averaged over the interval 1.45 < i^o < 1.47. In (a), rii = 50 and in (b), {Fg^l = 0.02. 
Full line: quantum result and dashed line: classical result. The arrows indicate the lifetimes of 
the ground state estimated by direct coupling to the continuum. The average is calculated on the 
logarithmic scale, i.e. shown is 10^^°Sio"^. 
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FIG. 5: Psuv averaged over the interval 1.45 < fo < 1-47. In (a), m = 50, in (b), {F^^'l = 0.02. 
Full line: quantum result, dashed line: classical result and dashed-dotted line: fit of classical result 
to a power-law. The arrows indicate the lifetimes of the ground state estimated by direct coupling 
to the continuum. The average is calculated on the logarithmic scale. 

We now compare the numerical signatures of quantum localization for intermediary times 
as seen in Psur, Figs. [T^ and 13 to the signatures seen in (n), Figs. |2t and El In all cases 
studied, the (n)^™ > (n)'^^ and a clear discrepancy between the quantum and the classical 
value can be seen for weak fields and low rij. The discrepancy decreases as I-Fq^I or rii 
increases, i.e. (n)''™ approaches the classical limit, however quite slowly. The approach to 
the classical limit appears faster in Pgur- Here signatures of quantum localization can only 
be see for weak fields and low rii while in the other cases shown, has already reached 
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the classical limit. We thus conclude that the more locally focused (n) provides stronger 
signature of quantum localization, i.e. the suppression of diffusion away from the initial 
state, than the Psur probing the entire bound part of phase space. 

B. Short-time dynamics: the cross-over to quantum localization 

We take now a closer look at the short time dynamics for the weak field strength. Fig. 
and Fig. |2b- Here a first surprise appears. One would, generally, expect PJjJj. and to 
agree with each other up to the localization time (or quantum break time, also referred 
to as Ehrenfest time) ti. The present case is non-generic in that the classical phase space 
distribution remains up to ti 200 kicks) more localized when one identifies Psur as a 
measure for localization. This quantum enhancement of ionization takes place even though 
the classical {n)ci moves closer to the ionization threshold while the quantum (n)qm remains 
close to the initial vale (Fig. |2b)- The origin is a true quantum effect: perturbative single 
photon absorption of high frequency itluq from higher-harmonic components with m > 
mc{ni) sufficient to reach the continuum. Here 

mcirii) = ni/(2z/o) . (11) 

A border for the field strength at which a single kick driving a high-lying Rydberg state 
displays quantum-classical correspondence can be found as follows j29| : The average classical 
energy transfer from a kick is 

AEk = Apl . (12) 

This energy transfer can be resolved quantum mechanically if AEk is not smaller than the 
quantum energy spacing AE'^. = l/rii, leading to the critical momentum transfer 

Ap?'* = (13) 

For Apo > Ap^^^, the quantum and classical distributions after a single kick agree all the 
way up to the threshold E = since the quantum level spacing decreases as n oo, 
and a good agreement between the quantum and classical survival probabilities is achieved. 
The difference between the classical and quantum distributions for Apo < Ap'^^^ can be 
understood by considering the dipole limit Apo < Apq'''°'^ = l/rii, for which the transition 
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operator for a single kick, exp{iAp q), reduces to 



exp {iAp q) ^ 1 + iApoq^Hi + ... . 



(14) 



In this limit, the quantum transition amplitude is proportional to Ap and the probability is 
proportional to Ap"^. By contrast, the classical probability is proportional to Ap^ implying 
that the quantum survival probability after one kick is smaller than the classical one j29|. 
Physically, this can be understood by considering the Fourier transform of a delta kick, 



indicating a "white" spectrum. The quantum system can absorb (virtual) photons of ar- 
bitrarily high frequency from the white spectrum accompanied by only a small momentum 
transfer Apo <^ 1. Processes with large energy transfer Ac^o ^ 1 but small momentum 
transfer Apo are far from the line Aujq ^ Ap'^/2 in the Ac^o — ApQ dispersion plane for 
a quasi-free electron, and thus effectively is inaccessible for a classical momentum transfer 
process. Only for deeply bound electrons with large local orbital momentum, po? near the 
nucleus, with Aujq ^ PqApq ^ Ap1/2 can such processes occur in the classical case. The 
density of classical phase space points for the initial state with such high po is, however, 
very small. In the quantum case, the high-frequency components in the driving field inter- 
act non-locally with the whole initial state leading to an enhanced ionization probability. 
This corresponds to dipole allowed transitions. Classical-quantum correspondence is only 
restored when classical diffusion in phase space dominates the quantum enhancement due 
to (virtual) photon absorption. 

For rii = 50 and uq = 1.45 the corresponding critical fields are |-Fo"*| = 0.033 and 
l^dipoioj _ Q gg47^ j_g_ ^Yie case studied in Figs. [T] and El are close to the dipole limit. The 
important point is that the dynamical role of the high-harmonics spectrum of (virtual) 
photons is responsible for the surprising, non-generic features of the periodically kicked 
Rydberg atom, different from other systems such as the Rydberg atom in a microwave field. 
This is directly verified by increasing the field such that the scaled momentum transfer 
ApQ ^ Apo"*, see Fig. IHl We show first a case in the transition region (I-Fq^I = 0.02, 
rii = 50) where still some traces of the non-generic features at short times can still be seen. 
For larger rii or larger l-Fg^^l such that Ap^ > Apg'''* the quantum-classical agreement for 
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FIG. 6: As Fig. Ht, however for (a): jF^'^j = 0.02 and = 50, (b): |F( 
and (c): \F^^\ = 0.02 and = 200. 
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short times is quite well fulfilled while discordance is found only for larger times, t > ti. The 
short-time dynamics now conforms with the naive expectation of close classical-quantum 
correspondence for all t < Ti. 

Well-known arguments j2fl| for the localization time lead to the order of magnitude esti- 
mate 

where (q) = is the typical spatial extension of the initial state, Aj = the wave length of 
the initial state, and (L) the mean Lyaponov exponent, ti is indicated by arrows in Fig. |2b 
and Fig. [3 For (n), ti gives a reasonable estimate for the number of kicks up to where 
quantum and classical dynamics mirror each other. 

The situation for the survival probability is more involved: For Apo < Apf^"*, Eq. ^ 
(see Figs. HI and Eti'), no agreement between P^^^ and is found for short times and a break 
time Ti now defined as the time where Pg^j. gets smaller than P^ due to quantum localization 
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Number of kicks 

FIG. 7: As Fig. Eb, however for (a): jF^'^j = 0.02 and rn = 50, (b): \F^^\ = 0.05 and rn = 50, 
and (c): \Fq'\ = 0.02 and = 200. The arrows indicate an approximate localization time defined 
in Eq. (fTT)|l. 

takes on a different meaning. Consequently, this break time is not well described by Eq. 

For the cases with Apo > ^Po^'*, (Figs. and Eb and c) we find no clear numerical signatures 

of localization in Psur- 

C. Long-time dynamics: the second cross-over 

Quantum localization in the kicked Rydberg atom is transient. After a large but finite 
number of kicks (Figs. ^ andlH)), the quantum survival probability rapidly decreases and 
at a second cross-over time (delocalization time) tj:, falls below even the classical value. 
Simultaneously, (n)qm drops sharply from values near rii to 1 (Figs.Eb andEb). Beyond the 
cross-over point, {n)c\ falls to values well below unity inaccessible to quantum mechanics. 
There, the residual fraction of the classical phase is "sheltered" and continues to decay slowly, 
i.e. algebraically. By contrast, the quantum bound-state probability decays exponentially in 
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the long-time limit. Beyond td the slow classical algebraic decay "wins" over the exponential 
decay. This novel scenario is markedly different from the sinusoidal driven Rydberg atom, 
where the quantum transport is bounded from below because of the regular phase space 
region, or, in a quantum picture, because of vanishingly small transition probabilities to 
states with low principal quantum numbers |31|. 

The rate of the long-time decay and, thus, of td can be estimated from the decay rate 
of the state n = 1 due to the high harmonics. This is the lowest-lying and most stable 
component of any coherent superposition forming a Floquet state. All frequencies of the 
driving field Eq. ^ with scaled energy muQ large enough to couple the Stark state having 
the largest overlap with the n = 1 hydrogenic state to the Stark continuum contribute to 
the transition probability according to Fermi's golden rule, 

= 2'Kp{E^)\F^''z{n = . (17) 

with Em = — 1/2 + muo/nf. Both dipole matrix elements z{n,E) and the density of con- 
tinuum states p{E) are numerically obtained by diagonalizing the Stark Hamiltonian (j21) 
in the pseudo-spectral basis. Summing over all contributions leads to a delocalization time 
To = (Xlm Tm)"^ which corresponds to within a factor of two to the lifetime of the most long- 
lived Floquet state with (n) ~ 1. {n)^^ averaged over an interval in u (1.45 < i^o < 1.47) 
together with the estimated td are shown in Fig. 01 td clearly gives a good estimate of 
the time-scale on which the break-down of quantum localization as seen in (n) takes place. 
The estimated td also coincide with getting smaller than the classical value (Fig. (HJ. 
We thus attribute the break-down of quantum-localization to the higher harmonics in the 
driving field. The breakdown of localization is expected if the ground state n = 1 is di- 
rectly coupled to the continuum. We note parenthetically that an experimental realization 
would require half-cycle pulses with high-frequency components in the UV region. Work is 
currently under way on a protocol to produce half-cycle pulses in the atto-second regime 

We now comment on the decay of the quantum survival probability for l-Fg^^l ~ 0.005 and 
0.02, Hi = 50, seen for intermediate times {K < 10^ and 10^ in Figs. and b, respectively). 
Here (n)'^"^ ^ rii (Fig. |3]) indicating the localization of the quantum distribution P{n,t)^^ 
(see also Fig. Ob)- As for (n), Eq. ((Zj), we attribute the decay of the locahzed bound-part 
of the quantum wavefunction in the kicked atom to the higher harmonics present in the 
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driving field, directly coupling the wave packet localized close to the initial state with the 
continuum. 



IV. FLUCTUATIONS IN THE SURVIVAL PROBABILITY 

The quantum survival probability shown in Fig. ^ display strong fluctuations under small 
variations of the kick frequency. Such fluctuations are a direct consequence of the photonic 



localization scenario 



1^, 



to be described in the following. 



A. High harmonics in the localization regime 

The observation of localization, i.e. the suppression of ionization, or, equivalently, the 
freezing out of portions of the wavefunction near n^, raises the question as to the underlying 
mechanism in the presence of the harmonic spectrum (Eq. )• High harmonics cannot 
only directly couple to the continuum, see Sect. 1111 CI but they also can lead to sequential 
excitation through a ladder of intermediate (quasi)bound states. This channel is the dom- 
inant mechanism for the excitation and ionization by the harmonic driving by the lowest 
harmonic uq. Jensen et al. [s] have discussed the suppression of the sequential ladder exci- 
tation by an harmonic driving as a mechanism for localization ("photonic localization"). It 
is therefore instructive to extend this approach to the present multi-photon case. 

Following 8] we assume quasi-resonant one-photon transitions to dominate the time evo- 
lution, i.e. we consider only transition between (bound) states with energy differences to 
the initial state approximately equal to k times the fundamental photon energy Eq = uo/rii. 
The detunings 

Ak = Ek- En, - ku (18) 

for the quasi-resonant states | k) form a pseudo-random sequence of numbers. Using the 
rotating wave approximation and setting Cfc(t) = bk{t) exp{—iAkt) , with bk{t) the expansion 
coefficient for the k*'^ quasi-resonant state in the interaction picture, leads to a set of coupled 
differential equations, 

'^ = J2HL'Ck' (19) 

k' 



with the matrix H'^ given by 



Hkk' = Skk' + Vk,k'- (20) 
14 



The semiclassical expression for the couphng matrix elements is j8| 

Vuk' ~ 0.411F/ ((nfcn,0'/'(mz/o)'/') (21) 

with m = \k — k'\ and Uk the main quantum number of the k*^ quasi-resonant state. While 
for the harmonically driven case Vkk' 7^ only for m = 1, the coupling matrix elements for 
the periodically kicked Rydberg atom in the m*^ off diagonal are proportional to m~^/^. Due 
to the randomness of the detunings A^, the matrix H'^ is thus a pseudo-random, power- 
law banded matrix. In random power-law banded matrices with the elements decreasing 
as m~^, the eigenstates are (weakly) localized for {3 > 1 [33]. Our extension of the model 
introduced in p] on the basis of [2] thus predicts that the eigenstates for the kicked Rydberg 
atom should be localized. Infering from the numerical observation of quantum localization 
that the detunings are "sufficiently" random, the generalization of the photonic localization 
scenario appears applicable. 

We finally comment on the (semi) classical border ^ 00. Increasing rii from 50 to 
200, the quantum localization for \Fq^\ = 0.02 almost vanishes (Figs. |3)d andEb)- This 
indicates that a delocalization border is present in the kicked atom similar to that found 
for the harmonically driven system [2j. The presence of a delocalization border implies a 
maximum for the applicability of the photonic localization theory. Detailed studies of 
this border in the kicked Rydberg atom remain to be performed. 



B. Characterization of parametric fluctuations 

We turn now to a quantitative description of the parametric fiuctuations of under 
variation of z/q seen in Fig. ^ Fig. |H1 displays the evolution towards an increasingly complex 
fiuctuation pattern as K increases. The amplitudes of the fiuctuations increase by several 
orders of magnitude. Strong fiuctuations in the survival probability have also been found 



e.g. in the harmonically driven Rydberg atom [ll| and in the kicked rotor when subjected 
to a varying Aharonov-Bohm fiux [34] or a variation of the kicking frequency . 

In order to characterize the increasingly finer scale on which these fiuctuations occur, we 
determine the average distance between adjacent maxima and minima on the frequency scale, 
{6uo{K)). At times smaller than the delocalization time td, {6uo{K)) is rapidly decreasing 
(Fig. EI). A similar build-up of fiuctuations in the survival probability with time is observed 
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FIG. 8: Survival-probability Pgur versus frequency v for some numbers K of kicks, = 50 and 



0.005. The dashed lines give the classical results. 



for the kicked rotor |35[. In the very- long time-limit beyond To, where the dynamics of 
the kicked system is governed by the decay of the ground state, no further fluctuations are 
produced but {6,y^^{K)) saturates. 
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FIG. 9: Average distance {6uq{K)) between extrema in Pgur as a function of number K of kicks. 
rii = 50 and |i^o^| = 0.005. We analyze a set containing 10^ points on 1.45 < z^o ^ 1.45008. 



It is now tempting to inquire whether the complexity of the fluctuations in the can be 
described by an approximate fractal dimension. For example, the fluctuations in the kicked 



rotor in the chaotic regime have been shown to have a fractal structure 134, hsl. Fractal 

n n 

conductance fluctuations have been predicted j36| and experimentally found, e.g. [37[, for 
the transport through cavities with mixed classical phase space. A fractal structure has also 
been found in the survival probability of a system with a mixed phase space 3^. In both 
cases, a semiclassical explanation based on on a classical power-law decay has been proposed. 
Since the kicked Rydberg atom displays classically such a power-law decay even though its 
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phase space is fully chaotic, rather than mixed, a fractal is a conceivable candidate. 
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FIG. 10: Fractal analysis of log]^o(-fsur) for n, = 50 and |-Fq^| = 0.005. In (a) we show the data 
analyzed for 2 x 10^ points in the interval 1.45 < z/q < 1.47 for different numbers K of kicks and in 
(b) we combine the wide set of (a) with data for a dense set containing 10^ points in the interval 
1.45 <uo< 1.45008 for K = 6.8 x lO'^ . 

For a given "resolution" Auq in frequency we calculate the local fractal dimension D(At'o) 
by means of a variational method j39|. Rigorously, this value should be independent of 
Auq. In the present physical context, we are satisfied if D(Az/o) is approximately constant 
over at least one order of magnitude in Az/q- The fractal analysis of the data shown in 
Fig. IHlis presented in Fig. ITUk. Since the survival probability fiuctuates over many orders of 
magnitude, we analyze the logarithm of the data. After more than 10^ kicks, a plateau starts 
to develop, the width of which reaches almost two orders of magnitude for larger K. The 
value of D on the plateau increases with K up to about D ^ 1.2. The process of increasingly 
finer rescaling is transient and stops beyond K > 10^ when T£, is reached. Zooming in on a 
narrow frequency interval for K > Kd, (Fig. ITOb). the plateau value is still is D ^ 1.2, but 
the width of the plateau is only about one order of magnitude. Thus, a tendency towards 
a non-integer dimension can be found in the positively kicked Rydberg atom. The value 
is weakly dependent on the number K of kicks for which the fractal analysis is made and 
appears to approach D ^ 1.2. 

It is now instructing to compare this value with the semiclassical prediction 
based on the power-law decay with exponent a = 1.5 (see Fig. ^ of the classical survival 
probability P^^^{K), Dsc = 2 — a/2 = 1.25. The value for Dsc is remarkably close to the 
plateau value found from the fractal dimension analysis. We thus conclude that the onset of 
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a self-similar fluctuation pattern can be observed with a dimension close to the semiclassical 
prediction derived from the classical power-law decay. We note that this process is transient 
in that the flnest scale for the fluctuations is determined by the time td beyond which the 
most stable bound state decays. 



V. SUMMARY 



We have studied the long-time limit of quantum localization of the positively kicked 
Rydberg atom, involving clear signatures of a quantum suppression of classical ionization. 
We compare the localization as seen in the survival probability to that seen in an aver- 
age quantum number (n) describing the position of the localized part of the wavefunction 
(see Eq. ((Tj) ). In (n) clearer signatures of quantum localization prevailing to higher fleld 
strengths and larger quantum numbers are found. Two cross-over times could be identifled. 
The cross-over from classical-quantum correspondence to localization (ti) and the destruc- 
tion or delocalization at a much later time {td)- Remarkably, beyond td, the quantum 
system decays faster then the classical counterpart due to direct transitions to the contin- 
uum resembling photoionization. This process, identifled in the present paper within a ID 
model, is expected to be operative in a full 3D model as well. Quantum localization is 
accompanied by strong fluctuations in the survival probability after a given number of kicks 
as a function of the frequency of the driving fleld. The average distance between the fluctu- 
ations {Suoit)) (X 1/t until t = Tf), whereafter (6^^) saturates. In the localization regime, the 
complex fluctuation pattern approaches a fractal. 

In this paper we have highlighted effects caused by the higher harmonics in the driving 
fleld distinguishing the kicked atom from the Rydberg atom driven by a microwave fleld. 
Further studies comparing these two systems, including an assessment on how closely the 
quantum localization in the kicked Rydberg atom is related to Anderson localization using 



e.g. the methods in j2 , 



is left for forthcoming studies. 
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